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Abstract 

We contrast four distinct versions of the BCS-Bose statistical crossover theory according to the form assumed for 
the electron-number equation that accompanies the BCS gap equation. The four versions correspond to explicitly 
accounting for two-hole- (2h) as well as two-electron- (2e) Cooper pairs (CPs), or both in equal proportions, or only 
either kind. This follows from a recent generalization of the Bose-Einstein condensation (GBEC) statistical theory 
that includes not boson-boson interactions but rather 2e- and also (without loss of generality) 2h-CPs interacting 
with unpaired electrons and holes in a single-band model that is easily converted into a two-band model. The GBEC 
theory is essentially an extension of the Friedberg-T.D. Lee 1989 BEC theory of superconductors that excludes 2h- 
CPs. It can thus recover, when the numbers of 2h- and 2e-CPs in both BE-condensed and noncondensed states are 
separately equal, the BCS gap equation for all temperatures and couplings as well as the zero-temperature BCS 
(rigorous-upper-bound) condensation energy for all couplings. But ignoring either 2h- or 2e-CPs it can do neither. 
In particular, only half the BCS condensation energy is obtained in the two crossover versions ignoring either kind of 
CPs. We show how critical temperatures T c from the original BCS-Bose crossover theory in 2D require unphysically 
large couplings for the Cooper/BCS model interaction to differ significantly from the T c s of ordinary BCS theory 
(where the number equation is substituted by the assumption that the chemical potential equals the Fermi energy) . 

PACS # 74.20.-z; 74.20.Mn; 05.30.Fk; 05.30.Jp 

Key words: Bose-Einstein condensation statistical model; BCS-Bose crossover theory. 

1 Introduction 

Boson-fermion (BF) statistical models of superconductivity (SC) as a Bose-Einstein condensation (BEC) [TJ [5] began 
to be seriously studied in the mid-1950's [3]-[6], pre-dating even the BCS-Bogoliubov statistical theory [7]-[9]. Although 
BCS theory only contemplates the presence of "Cooper correlations" of single-particle states, BF models [3]-[6] [1"0]-[T8] 
posit the existence of actual bosonic CPs. A drawback of early BF models is the notorious absence of an electronic gap 
A(T), with T the absolute temperature. Perhaps the first BF model with a gap was introduced in Ref. [TH]- Somewhat 
later, the remarkable relation A(T) oc y/no(T), with no(T) the BEC condensate number density of electron pairs, first 
seems to have appeared [20 . It resurfaced a year later in the BEC BF model in 3D of Friedberg and T.D. Lee [HI |2"2"] 
applied to cuprate superconductors. With just one adjustable parameter (the ratio of perpendicular to CuCVphine 
boson masses) this theory fitted (22] quasi-2D cuprate T c /Tp empirical values [23] rather well. The ratio turned out 
to be 66,560 — just under the 10 5 anisotropy ratio reported for B2+ x Sr2- y Cu02±s [24] almost contemporaneously. 



An extension of the work in Refs. [2TJ [22] is a generalized BEC (GBEC) statistical single-band theory whereby 
a superconducting BCS condensate was recently suggested [10] , and subsequently confirmed (but only to lowest 
order in the BF coupling) to be precisely a Bose-Einstein condensate (BEC) of equal numbers of bosonic two-electron 
(2e) and two- hole (2h) Cooper pairs (CPs), at least inasmuch as the GBEC reproduces the same BCS gap equation for 
all temperature and coupling as well as the same T = condensation energy found from BCS theory. The holes that 
make up the 2h-CPs originate precisely the Fermi sea asssociated with the iV-electron system in the simple single-band 
model studied here. One advantage of the single-band model is that it allows recovering, among other theories, the 
BCS theory as a special case. The distinction (Ref. [35] pp. 70-72) between single particles and holes is, in a sense, 
trivial. Not so for particle-pairs and hole-pairs, as will be seen shortly. 

The BF coupling assumed appears in an interaction many-body Hamiltonian Hi nt which defines the GBEC theory. 
Added to Hi nt is an unperturbed Hamiltonian Hq describing a free ternary gas of unpaired electrons plus 2e-CPs plus 
2h-CPs. The noninteracting ternary gas represents the normal state of the original, strongly-correlated many-electron 
system under study, and is a viable candidate for a so-called "non-Fermi-liquid." The new GBEC theory embodied in 
H = Hq + Hi n t is in essence a complete BF (statistical) single-band model that, however, admits departure [11] from 
the perfect 2e-/2h-CP symmetry that constrains BCS theory by construction. It can be diagonalized via a Bogoliubov 
canonical transformation exactly if one neglects nonzero center-of-mass- momentum (CMM) CPs in Hi nt as is done in 
BCS theory but not in Ho which in BCS theory represents a pure electron gas. The GBEC theory is complete only 
in the sense that 2h-CPs are not ignored. It reduces to all the known statistical theories of superconductors (SCs), 
including the BCS-Bose "crossover" picture in the four versions to be distinguished below. Its practical impact is 
that it yields [121 1 1 3] robustly higher T c 's than BCS theory without abandoning electron-phonon dynamics when one 
departs from the perfect 50/50 symmetry of 2e-/2h-CPs in the condensate. 

In the literature, electron-phonon dynamics have been widely mimicked by the s-wave BCS/Cooper model interac- 
tion Vk^' It is a nonzero negative constant — V, if and only if single-particle energies lie within the energy 
interval [max{0, (i — huo}, /i+Swd] where fj, is the electron chemical potential and ujd the Debye frequency. We em- 
ploy this model interaction here. Other pairing symmetries beyond pure s-wave can also be accommodated. Although 
it sheds considerable light on different possible crossover schemes, our single-band picture where single-electron and 
single-holes are assumed to have the same effective mass is not as realistic in describing real materials as a multiband 
(say, a valence-like band for holes and a conduction-like one for electrons) theory where these two masses can differ, 
as will be discussed later. 

A fundamental drawback of early [3]- [6] BF models, which took 2e CPs as analogous to diatomic molecules in a 
classical atom-molecule binary gas mixture, is the cumbersome introduction of an electron energy gap A(T). "Gapless" 
models, however, are useful [151116] in locating transition-temperature singularities if approached from above, i.e., from 
the normal state where T > T c . 

The "ordinary" CP problem [26] for two distinct interfermion interactions (the (^-potential well [TTl [18] or the 
Cooper/BCS model [7] [26] interactions) neglects the effect of 2h-CPs treated on an equal footing with 2e [or, in 
general, two-particle (2p)] CPs. On the other hand, Green's functions [25] can naturally deal with hole propagation 
and thus accommodate both 2e- and 2h-CPs via, e.g., the Bethe-Salpeter equation (27] [28]. In addition to the 
generalized CP problem, a crucial result [101 E] as already mentioned is that the BCS condensate consists of equal 
numbers of 2e- and 2h-CPs. This was implicitly already suggested from the perfect symmetry about electron energy 
e = /i of the well-known Bogoliubov [29] v 2 (e) and u 2 (e) coefficients, with the tail of v 2 (e) above e — n representing 2e 
correlations and that of u 2 (e) below e — fi refers to 2h correlations. 

In this paper we show how: a) four versions of the BCS-Bose statistical crossover theory can be obtained by ignoring 
either 2h- or 2e-CPs or by including both; c) for only two of the four versions can the precise BCS gap equation for 
all temperatures T be derived; c) crossover-picture T c s, defined self-consistently by both the gap and fermion-number 
equations, requires unphysically large couplings (at least for the Cooper/BCS model interaction in 2D SCs) to differ 
significantly from the T c of ordinary BCS theory defined without the number equation since here the chemical potential 
is assumed equal to the Fermi energy; and d) the full T = BCS condensation energy follows from one crossover 
version but only half of it from the two versions ignoring either kind of CPs. The condensation energy is simply 
related to the ground-state energy of the many-fermion system, which in the case of BCS is a rigorous upper bound 
to the exact many-body value for the given Hamiltonian as BCS theory starts from a variational wave function for 
the superconductor ground state. These results, with the exception of (c) which does not apply, are expected to hold 
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also for neutral- fcrmion supcrfluids (SFs) — such as liquid 3 Hc [30, 31 , neutron matter and trapped ultra-cold fermion 
atomic gases |32j-|43) — where the pair-forming two-fermion interaction, of course, differs from the Cooper/BCS one 
for SCs. 

2 Generalized BEC Theory (GBEC) 

The GBEC theory is described in detail in Refs. |10j-[13j: here we summarize its main equations. It applies in d 
dimensions and is defined by a Hamiltonian of the form H = Hq + Hi nt . The unperturbed Hamiltonian Ho should 
ideally be, to quote Leggett [44] "an appropriate 'zeroth-oder' starting point" accounting for "pairs of electronic 
excitations with charge 2e that all have the same ground-state wavefunction." Thus, our Ho corresponds to a non- Fermi- 
liquid "normal" state which, besides just fermions, is an ideal (i.e., noninteracting) ternary gas mixture of unpaired 
fermions and both types of CPs namely, 2e and 2h, the latter introduced without loss of generality. Specifically 

H = J2 ^a+^at^ +Y,E+(K)b+b K -Y,E-(K)c+CK (1) 

ki.si K K 

where K = ki + k 2 is the CMM wavevector of the pair, while = ti 2 kf/2m are the single-electron, and E±(K) 
the 2e-/2h-CP phenomenological, energies. Here a£ (flki.si) are creation (annihilation) operators for fermions and 
similarly 6 K (6k) and c K (ck) for 2e- and 2h-CP bosons, respectively. These b and c operators depend only on K and 
so are distinct from the BCS operators depending on both K and the relative k = ^(kj — k 2 ) discussed in Ref. [7] 
Eqs. (2.9) to (2.13) for the particular case of K = and shown there not to satisfy the ordinary Bose commutation 
relations. But because two pairs cannot exactly overlap in real space without violating the Pauli Principle, they are 
often considered "hard-core bosons," albeit of hard-core radii + . For this reason, one can probably not expect to 
be able to construct the b and c operators directly from the a operators in order to establish that b and c obey Bose 
commutation relations precisely. Nonetheless, these pairs stand for objects that can easily be seen to obey Bosc- 
Einstcin statistics as, in the thermodynamic limit, an indefinitely large number of k values correspond to a given 
K value defining an energy level E + (K) or E_(K), This is all that is needed to ensure the BEC (or macroscopic 
occupation of a given state that appears below a certain fixed T = T c ) found numerically a posteriori in the 

GBEC theory. Furthermore, being non- interacting (except for the Pauli Principle restriction mentioned), CPs satisfy 
the Ehrenfest-Oppenheimer [45] criteria for two clusters of charges to conserve a specific kind of statistics, either Bose 
or Fermi. These assumed properties are justified a posteriori when in the GBEC theory: a) the BCS gap equation is 
recovered for equal numbers of both kinds of pairs, both in the K = state and in all K ^ states taken collectively, 
and in weak coupling, regardless of CP overlaps; and b) the precise familiar BEC T c formula emerges [TU] when i) 
2h-CPs are ignored, the Friedberg-T.D. Lee model [21] [22] equations are recovered and ii) one switches off the BF 
interaction. The only difference in the recovered BEC T c formula is that the boson number density now depends on T c , 
as expected in a boson- fermion mixture where populations are T-dependent. Finally, we note that fermion scattering 
terms [46] are not included in [T] as they are not expected to be substantial, say, in the BCS limit of high electron 
density where they would be the most effective, which in turn is included in the GBEC model as a special case. 

Two-hole CPs in ([T]) are postulated to be distinct and kinematically independent from both the 2e-CPs and the 
unpaired electrons, i.e., operators a, b and c are assumed to commute with each other. This postulate is grounded on 
magnetic- flux-quantization measurements establishing the presence of pair charge carriers in both conventional [47[ 148] 
as well as cuprate [49] superconductors, and on the fact that no experiment has yet been done, to our knowledge [50j . 
that distinguishes between electron and hole CPs. The latter uncertainty further motivates a Hamiltonian such as ([T]) 
with both kinds of CPs. 

The interaction Hamiltonian Hi nt in the expression H = Ho + Hi nt describes the formation and disintegration of 
CPs, respectively, from and into unpaired electrons and holes. It is further simplified by dropping all K ^ terms. 
This is also done in BCS theory in its full Hamiltonian H — Hq + Hi nt , but kept in the GBEC theory in its unperturbed 
Ho portion JT]). The GBEC Hi nt is made up of four distinct BF interaction vertices each with two-fermion/one-boson 
creation and/or annihilation operators. These vertices depict how unpaired electrons (subindex +) [or holes (subindex 
— )] are involved in the formation and disintegration of the 2e- (and 2h-) K = CPs in the d-dimensional system of 
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size L, namely 



H mt = L- d ' 2 f+(k){a+ r a± k L b + a_ ka a k , T &+} + ZT d/2 ^ /-(fc){<, T a+ ka 4 + 0-k,iO klT Co} (2) 

k k 

where k = ^ (k x — k2) is again the relative wavevector of a CP. The interaction vertex form factors f±(k) in ^ arc 
essentially the Fourier transforms of the 2e- and 2h-CP intrinsic wavefunctions, respectively, in the relative coordinate 
of the two fermions. The GBEC theory is thus reminiscent of the Sommerfeld theory of the electron gas combined 
with the Debye picture of the phonon gas which together give a binary mixture of noninteracting electrons and 
phonons, a picture which describes low-T specific heats in metals and insulators. But to explain either resistance and 
superconductivity, they must then be allowed to interact via the Frohlich electron-phonon interaction [51j of a form 
analogous to ^ but without hole terms. In contrast, the full BCS Hamiltonian Hq CS + Hf^ s consists of only the 
first (electron) term on the rhs of ((T|), namely 



- ^kiKus^^i,^ (3) 



tjBCS 

ki,si 



and 

H?r?t S = ^A^-^-hiaht (4) 

ki,li 

The BCS Hq CS thus represents a Fermi liquid normal state. 

The interaction vertex form factors f±(k) in are essentially the Fourier transforms of the 2e- and 2h-CP intrinsic 
wavefunctions, respectively, in the relative coordinate between the paired fermions of the CP. In order to eventually 
recover BCS theory, in Refs. [TD]- [T3] the corresponding energy form factors were picked as 

, ._( f for E f < e < E f + Se , . 

J + { ' \ otherwise, { ' 

|/ for E f -5e<e<E f 
w [0 otherwise. v ; 

This is after one introduces the quantities Ef and Se as new phenomcnological dynamical energy parameters (in 
addition to the positive BF vertex coupling parameter /) that replace the previous phenomenological CP energy 
parameters E±(0), through the definitions 

E f = ±[#+(0) + £L(0)] and Se = ±[E+(0) - £L(0)] > (7) 

where -E+(0) and E-(0) are the (empirically wnknown) zero-CMM energies of the 2e- and 2h-CPs, respectively. Note 
that 2Ef lies midway between E + (0) and i5_(0). Alternately, instead of ([7]) one can write the two relations 

E±(0) =2E f ±5e. (8) 

The quantity Ef serves as a convenient energy scale; it is not to be confused with the Fermi energy Ep = ^mvp = kgTp 
where Tp is the Fermi temperature. The Fermi energy Ep equals Trh 2 n/m in 2D and (S 2 /2m)(37r 2 n) 2 / 3 in 3D, with 
n = N/L d the total number-density of charge-carrier mobile electrons, while Ef is of the same form but with n 
replaced by, say, rif, which in turn serves as convenient electron-density scale. The quantities Ef and Ep coincide 
only when perfect 2e/2h-CP symmetry holds, i.e., when n = nf. 

The grand potential f2 for the full Hamiltonian H = Hq + Hi nt given by |T]) and ((2|) is then constructed via (Ref. 
[25] Eq. 4.14) the definition 

n(T,L d ,fi,N ,M ) =-k B T\nWve-^ H -^A (9) 

where "Tr" stands for "trace" and (3 = 1/ksT with T the absolute temperature. It is related to the system pressure 
P, internal energy E and entropy S by — —PL d = F — fiN = E — TS — fiN, where F is the Helmholtz free energy. 
Following the Bogoliubov prescription [52] , one sets , bo equal to y/No and Cq , cq equal to \/A1q in , where No is 
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the T-dependent number of zero-CMM 2e-CPs and Mo likewise for 2h-CPs. This allows exact diagonalization for any 
coupling, through a Bogoliubov transformation of the a + , a fermion operators, giving [53] after some algebra 

/>oo />OC 

n(T,L d ,[i,N ,M )/L d = / deN{e)[e- n-E(e)}-2k B T / deN{e) ln{l + exp[-/3E(e)]} 

Jo Jo 

/>oo 

+ [E+{0) - 2fi]n + k B T / deM(e)ln{l-exp[-(3{E + (Q)+e-2^}}} 

Jo+ 
1*00 

+ [2/j,- E-{0)]m o + k B T / deM (e) ln{l - exp[- /3{2/i - E_(0) + e}}}. (10) 
Here N(e) and M(e) are respectively the electronic and bosonic density of states, while 



E(e) ee V(e-M) 2 + A2( £ ) = ^/ (e - ^ + n fl(e) + m p_(e) (11) 

since A(e) = ^/mf + {e) + ,/m^/_(e) and /+(e)/+(e) = from © and ©, with n (T) = N (T)/L d and m (T) = 
Mo(T)/L d being the 2e-CP and 2h-CP number densities, respectively, of BE-condensed (i.e., with K = 0) bosons. 

Minimizing F with respect to Nq and Mq, while simultaneously fixing the total number N of electrons by intro- 
ducing the electron chemical potential /1 in the usual way, namely 

9F 0, ™ = 0, and ^ = -AT (12) 



cW ' dM Oil 

ensures an equilibrium thermodynamic state of the system with volume L d at temperature T and chemical potential 
fi. Evidently, N includes both paired and unpaired CP electrons. Some algebra then leads [53] to the three coupled 
integral Eqs. (7)-(9) of Ref. [TO] which, since from ([3]) and ([5]) one has that / + (e)/_(e) = 0, can be simplified to the 

two "gap-like equations" 

Ef+5e 



1, 2 f ^ ATf ^ tanh \Py/[e- /x(T)] 2 + f 2 n (T) 



p% + fe-2MT)]=2/ J / ^iV(e) ^ ^ - (13) 

2 J y/[e- fi{T)} 2 + f 2 n {T) 



E 



1 a f ^ Ar/ ^tanh - n(T)] 2 + f 2 m {T) 



[2^T)-2E f +Ss] = -f / deN(e) ^ " " L=p (14) 

Ef—5e 

and a single "number equation" (that guarantees charge conservation) 

2n s (T)-2m B (T)+n, / (T) = n. (15) 

where 

n/(T) = | deiV(e)[l - ^ tanh l/325(e)] (16) 


is clearly the number of unpaired electrons. This is identical with 2^ k v 2 ,(T) with v 2 (T) the well-known T-dependent 
Bogoliubov -^-coefficient. In (fT5"]) n = N/L d is the number density of electrons while n B (T) and m B (T) are, respec- 
tively, the number densities of 2e- and 2h-CPs in all bosonic states (both K = as well as K > 0). The "complete" 
number equation (|15p can be rewritten more explicitly as 

2n Q (T) + 2n B+ (T) - 2m (T) - 2m B+ (T) + n f (T) = n crossover version A (17) 

where n B (T) is [13] 

00 

n s (T) = n (T) + n B+ (T); n B +(T) = / cfeAf ( £ )[exp/3{£; + (0) + e - 2 A t} - lp 1 (18) 

0+ 
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and similarly for mg(T) which is 



m B (T) = m (T) + m B+ (T); m B+ (T) = J deM(e)[exp /3{2fi-E_(Q) + e} - l]" 1 . (19) 

d+ 

Clearly, m B+ (T) are precisely the number of "pre-formed" K > 2h-CPs, and n B +(T) that of 2e-CPs. These CPs are 
non-condensed in contrast with the K = CPs which are BE condensed. Evaluating the integrals requires knowing 
the bosonic density-of-states M(e) of CPs of energy e, which in turn requires knowing the dispersion relation e vs. K, 
e.g., as has been determined via the Bethe-Salpeter equation in the ladder approximation in 3D [27) and in 2D [28j . 

Self-consistent (at worst, numerical) solution of the three coupled integral equations (flU)) . (|14j) and (jTTJ) then yields 
the three thermodynamic variables of the GBEC theory 

no(T,n, /i), mo(T,n, fx), and (x(T,ri). (20) 

Figure 1 displays the three BE condensed phases — labeled s+, s— and ss — along with the normal phase n, that 
emerge [TT] from the GBEC theory. Phase s+ stands for a pure 2e-CP BE condensate, s— for a pure 2h-CP such 
condensate and ss denotes a mixed phase. Only the two pure phases were found [11] to display T c values higher than 
the corresponding BCS value, while the mixed phase occurs below this value. 




Figure 1: Illustration in the no-mo plane of three GBEC theory condensed phases (the pure 2e-CP s+ and pure 
2h-CP s— BE condensate phases and a mixed phase ss) along with the normal (ternary BF non-Fermi-liquid) phase 
n that corresponds to the origin at no = = too- 

For the two pure phases one can, in principle, shift from the single-band model implied so far to a two-band model 
by allowing the particle (e) masses to differ from hole (h) masses; this can be done by introducing two different Fermi 
energies E e F and Ep that differ precisely by these two masses. 

The GBEC theory contains [13] the key equations of all five distinct statistical theories as special cases. These 
range from ordinary BCS to ordinary BEC theories, which are thereby completely unified by the GBEC theory. 
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Perfect 2e/2h CP symmetry signifies equal numbers of 2e- and 2h-CPs, more specifically, ns(T) = ms(T) as well as 
no(T) = 7tiq(T). This implies that Hb+(T) = mB+{T) for all T, meaning that the exponents in (|18|) and (|19p coincide 
so that with |5]| this makes Ef = \i. The GBEC theory then reduces to the gap and number equations [viz., in 2D 
for T = T c both (|28[) and (f2"5)) below] of the original [55] BCS-Bose crossover picture with the Cooper/BCS model 
interaction — if its parameters V and Hlod are identified with the BF interaction GBEC Hamiltonian Hi nt parameters 
f 2 /2Se and Se, respectively. This one-to-one correspondence between Hi nt and H^ s defined in ([2]) and ([4} justifies 
the particular choice of form factors ([5J and © for the BF interaction. The original crossover picture for unknowns 
A(T) and /x(T) is now supplemented by the central relation 

A(T) = fy/no(T) = fy/mo(T). (21) 

All three functions A(T), uq(T) and mo(T) have the familiar "half-bell-shaped" forms. Namely, they are zero above 
a certain critical temperature T c , and rise monotonically upon cooling (lowering T) to maximum values A(0), no(0) 
and mo(0) at T = 0. The energy gap A(T) is the order parameter describing the superconducting (or superfluid) 
condensed state, while hq(T) and mo(T) are the BEC order parameters depicting the macroscopic occupation that 
occurs below T c in a BE condensate. This A(T) is precisely the BCS energy gap if the boson- fermion coupling / is 
made to correspond to \J2VTujjd within the GBEC formalism. Evidently, the BCS and BEC T c s are the same. Writing 
(I2ip for T = and dividing this into (|2ip gives the much simpler /-independent relation involving order parameters 
normalized to unity in the interval [0, 1] 

A(T)/A(0) - y/MT)/M0) = Vm o (T)/m o (0) — ♦ 1 

0. (22) 

T>T C 

The first equality, apparently first obtained in Ref. [20] . connects in a simple way the two heretofore unrelated 
"half-bell-shaped" order parameters of the BCS and the BEC theories. The second equality [TUl [11] implies that a 
BCS condensate is precisely a BE condensate of equal numbers of 2e- and 2h-CPs. Since (|22|) is independent of the 
particular two-fermion dynamics of the problem, it can be expected to hold for either SCs and SFs. 



3 Gap equation 

The standard procedure in all SC and SF theories of many-fermion systems is to ignore dealing explicitly with 2h-CPs 
altogether. Neglecting in (110|) all terms containing mo(T), E_(0) and /-(e) leaves an Q(T, L d , fi, No) defining a binary, 
instead of ternary, BF model. Minimizing the associated Helmholtz free energy F(T, L d ,/d,No) = 0(T, L d , //, No) + fiN 
over iVo (for fixed total electron number N) requires that dF/dNo = = dF/dno, which becomes 



f 

Jo 



deN(e) 



2exp{-/3£;(e)} 



l + cxp{-/3 J B(e)} 



fM + [S +( 0)-2 M ] = 
dno 



or 



2[2E f + 5e-2[i} = f 2 deN{e)— —tarih- (3 E{e). (23) 

JE f E \ e ) 2 

Using ijHJ) yields precisely the BCS gap equation for all T, Eq. (3.27) of Ref. [7], provided one picks Ef — fi, namely 



1 = A / d£ - = tanh -[3Je + A 2 (T) (24) 

J V^ 2 + A2 ( T ) 2 

where ^ = e — /x, since A = N(Ep)V = f 2 N(Ep)/2Se while fe = [see relation between V and / stated just above 
(|21[)]. and provided iV(e) can be taken outside the integral sign in (f23|) . This last operation is exact in 2D when N(e) 
is independent of e and is otherwise a good approximation if Hu>d <C /i. 

However, the choice Et = [i cannot be justified, to our knowledge, without assuming within the GBEC that 
ub(T) = uib(T) as well as uq(T) = mo(T), i.e., by explicitly recognizing the existence of 2h-CPs along with 2e-CPs 
and taking them in equal or 50-50 proportions. 
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4 Number equation 



Besides the normal phase consisting of the ideal BF ternary gas described by Hq, three different stable BEC phases 
emerge [11] when solving all three equations (fl~3]l to (fl4|) or (fT7|) : two pure phases, a pure 2e-CP BEC and a pure 
2h-CP BEC, as well as a mixed phase consisting of both types of BECs in varying proportions. For a half-and-half 
mixed phase, i.e., i%o(T) = mo(T) and ub+(T) = rriB+(T), all the boson number-density terms in (|17p cancel and the 
BCS number equation 

n — nj(T) crossover version B (special case of A) (25) 

is recovered, with n/(T) defined by (|16p . Crossover version B does not explicitly neglect either kind of CP, nor does 
it draw a distinction between them, and is the version applied below in 2D to obtain (|2"8"|) to ([50)) and the results of 
Fig. 2. This is the original crossover version first presented in 1967 in Ref. [55] , 

If 2h-CPs are ignored altogether, the companion number equation follows from the last equation of (fT2]) as 

n = rif(T) + 2ub{T) crossover version C (26) 

where rif(T) is interpreted as the number density of unpaired but BCS-correlated electrons and is given by (|16[) . In 
perhaps the first attempt to discuss [TI5] BEC in 1969 within the BCS-Bose crossover picture, Eagles [S3] imposed 
(|26p to accompany the gap equation in what was perhaps the first BF model with a gap. This differs from the much 
simpler number equation of crossover version B, which gave ([2"§1) below as a special case for T = T C when A(T C ) = is 
substituted into ([25]) and (fl6]) if one uses the identity 1 — tanh(ir/2) = 2/ (exp x + 1). It is this version that corresponds 
to the Friedberg-T.D. Lee model [2T]l22]. 

Similarly, ignoring 2e-CPs and keeping only 2h-CPs leads to f2(T, L d , /i, Mo) from which to minimize F(T, L d , /i, Mo) 
over M requires that one set dF/dMo = = dF/dm . Noting that £*(£) = E(—£), this also leads to the gap equation 
(1241) provided, again, one picks Et = fi, but now with the companion number equation 

n = n/(T) — 2tob(T) crossover version D (27) 

instead of ([2"6"]) but with the same n,f(T) as in (TT6"]) . 



5 BCS-Bose crossover T c compared with BCS T c in 2D 

The original crossover theory [55j is defined by two simultaneous coupled equations, the BCS gap and number equations, 
without the BCS assumption that the chemical potential /i equals the Fermi energy Ep. For subsequent extensions of 
the original version, see reviews in Refs. J56][57]- The critical temperature T c is defined by A(T C ) = 0, and is to be 
determined self-consistently with /i(T c ) by solving both gap and number equations. Because of its interest in quasi-2D 
cuprate superconductors [55], in this section we concentrate on 2D only. For the Cooper/BCS model interaction, if 
A = N(Ep)V where N(Ep) = m/2Trh 2 , the two crossover equations to be solved self-consistently reduce to 

,ftu, D /2fc B T c tanhx rhu D /2k B T c tanha , 
1 = A / dx (if [i > huj D ); 1 = A / dx— (if /.i < hu D ) (28) 



M(T c )/2fc B T c 



2.r 



' U = E F . (29) 



/o exp{[e-M(7 1 c)]/2fci 3 T c } + l 
The last integral can be done analytically and leaves the explicit expression 

M (T c )=A; B T c ln( e ^/feBT c _ 1) (30) 

which is then eliminated symbolically from ([28]) to render T c as an implicit function of A alone. Using hujp>/Ep = 0.05 
as a typical value for cuprates, increasing A makes /J.(T C ) decrease from its weak-coupling (where T c — > 0) value of Ep 
down to hwo when A ~ 56, an unphysically large value as it well exceeds the Migdal ionic-lattice stability upper limit 
[59] of 1/2, although [60] "there is no universally accepted, simple, and quantitative stability criterion." 

Fig. 2 displays T c (in units of Tp) as function of A. Note that room-temperature SCs (RTSC) are predicted by 
BCS-Bose crossover theory but only for A values definitely larger than about 10 that are still too unphysical. 
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Figure 2: Critical SC temperatures T c in units of Tp for the BCS-Bose crossover theory (full curve) in 2D compared with 
the BCS value from the exact implicit T c equation (see, e.g., Ref. [25], p. 447) 1 = A J^ UJD ' 2kBTc d xx ~ l tanhrr (upper 
dashed curve) valid for any d > in any coupling A, and its weak-coupling explicit solution T c ~ 1.134/iwd exp(— 1/A) 
(lower dashed curve). The dot-dashed "appendage" signals a breakdown in the BCS/Cooper interaction model when 
n(T c ) turns negative, as the Fermi surface at fi then washes out and the interaction model becomes meaningless. The 
value of A = 1/2 marked is the maximum possible value allowed just short of lattice instability in 3D for this interaction 
model, at least by one criterion [55] ■ Unshaded region refers to room-temperature superconductivity (RTSC) for SCs 
with Tp < 10 3 K. 



6 Condensation energy 

The T = condensation energy per unit volume according to the GBEC theory, given QlUp. is 

e s - En = n fl (r = o)-n n (r = o) 

L d ~ L d 



(31) 



since for any T the Helmholtz free energy F = E — TS = fi + ^N, with S the entropy, and n is the same for either 
superconducting s or normal n phases with internal energies E s and E n , respectively. In the normal phase uq(T) = 0, 
m Q (T) = so that A(T) = for all T > 0, so that (TO]) reduces to 

" n( ^T 0) = J™ deN(e)(e - M - \e - »\) = 2 £ deN(e)(e - M ) = 2 J_ d£N(QZ. (32) 

For the superconducting phase, and when no(T) — mo(T) and ub{T) — rns(T) hold, i.e., crossover scenario B, one 
deduces from (SJ and (fTUj) that /j, = Ef. Putting A(T = 0) = A in ([T0| as well as Se = fnoo, while using ([8]), gives 



^ d 0) = 2^n (0) + |°° ^iV(C) (€ - V^+A^) 
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= 2hw D n (0) + 2 



—hur 



tUOr 



(33) 



The first factor of 2 in the last line comes precisely from the condition tiq(T) = mo(T) while the last two factors of 2 
arise from the condition that according to ([5]) and ([6]) the magnitudes of /+(e) and /-(e) are the same and equal /. 
Subtracting ([32]) from ([33]) and putting iV(£) = -ZV(O), the density of electronic states at the Fermi surface [designated 
before as N(E F )] yields 

77 rp rTuiJD , 

" : ' --2huj D n {0) + 2N(0) dH£ 



2huj D n {0) +N(0) 



L d 



A 2 In- 



ve + A 2 

A 



(GBEC) (34) 

V(^d) 2 + A 2 _ 

for T = gives Eq. (2.40) of Ref. [7], 

(35) 



exactly, by standard integrations. Using the expression that follows from 
namely 

_ huj D 
~ sinh(l/A) 

where A is related to GBEC BF interaction parameter / through 

A = VN(0) = f 2 N(0)/2hio D . 

This makes the first term on the rhs of (f3~4"]) exactly equal to A 2 7V(0)/A which in turn can be shown to cancel exactly 
against the log term if one recalls the hyperbolic-function identity sinh 2 x + 1 = cosh 2 x. Thus, the GBEC theory 
condensation energy (fM]) is identical for any coupling to that of BCS theory, Eq. (2.42) of Ref. [7J, namely 



E s — E n 
L d 



N(0)(Hljd) 



1- A/1 



-7V(0)A 2 



1 

1 - - 

4 



{A/huj D 
A 



D 



O 



(BCS) 
A 



(36) 



ftwi 



This energy, associated with the expectation value of the BCS trial wavefunction gives a rigorous upper bound to 
the exact ground-state energy of the BCS Hamiltonian. Empirically, for niobium (Nb, bec, T c ~ 9.3-ftT, critical 
magnetic field H c ~ 160k A/ m) the condensation energy to be compared with the BCS result ([36]) works out to be 
just 2 x 10 _6 eV/atom [62]. The equivalence of (l34|) and (|36|) seems to suggest that, as in the GBEC theory, there are 
no pair-pair interactions in the BCS theory either, as is evident from Hamiltonians |T]), ([2|) and Q. 

What happens on ignoring either 2e- or 2h-CPs, as seems to be common practice in theories of SCs and SFs? This 
gives crossover versions C and D. Starting from (fTU)) for T = 0, and following a similar procedure to arrive at (f3"3"]) but 
without 2h-CPs such that /_ = 0, to o (0) = and n (0) = A 2 // 2 , one gets 



n s {T = Q) 



Subtracting 



L' 1 

from ([57)) gives 
E. 



= huj D n (0) 



d£N(t)t + N(0) 



Kit- ve+A* 



E„ 



fuu D n o (0) 



+ n(o) di (z - veT^) 



(37) 



(38) 



which is just half the full GBEC theory result (f3"4"]l . Furthermore, if [(E s — E n )/L d ]- is the contribution from 2h- 
CPs alone we may assume that / + = and tiq(O) = and eventually arrive at precisely the rhs of ([38]) but with 
m (0) = A 2 // 2 in place of n (0) = A 2 // 2 . Hence 



E s — E n 




E s — E n 


L d 


+ 


L d 



1. - -iV(0)A^ 

a^o 4 



1 / A 



4 \ fkoi 



1- \ 1 



O 



f {A/Hujd) 

(— 

\huj D 



(39) 
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which again is just one half the full GBEC theory condensation energy (j34|) that was found to be identical to the 
full BCS condensation energy (J3SJ) - Though not too surprising as the function E(e) = \/(e — /i) 2 + A 2 (e), where 
A(e) = ^/n~of+(e) + v / mo/_(e), becomes "half-gapless" in either crossover versions C or D, this one-half difference 
occurs precisely because either no or mo, and /+(e) or /-(e), have been deleted. Including both 2e- and 2h-CPs 
gave similarly striking conclusions on generalizing via the Bethe-Salpeter equation the ordinary [26] CP problem from 
unrealistic infinite- lifetime pairs to the physically expected finite-lifetime ones of Refs. [27l [28] . 

7 Conclusions 

The recent generalized BEC (GBEC) statistical single-band theory was employed to distinguish four different versions 
of the BCS-Bose crossover picture. One of these is the original BCS-Bose crossover theory with number equation 
([25]) . crossover version B. For the Cooper/BCS model interaction predicts in 2D virtually the same T c s to well beyond 
physically unreasonable values of coupling, as the (allegedly less general) BCS statistical theory where the number 
equation becomes trivial on assuming that the electron chemical potential u — Ep, the Fermi energy. However, T c s 
much higher than those of the BCS-Bose crossover theory have been obtained [53] via the GBEC number equation 
(TTT]) . designated here as crossover version A, that includes both electron- or hole-pair bosons explicitly but in different 
proportions. 

The GBEC statistical theory also reveals that the BCS gap equation for all temperatures follows rigorously only 
when neither hole- nor electron-pairs are ignored and occur in equal proportions, separately for zero- and nonzero- 
CMM pairs, and that the resulting GBEC T = condensation energy equals the entire (rigorous-upper-bound) BCS 
value for any coupling. But that it is only half as large when either kind of pair is ignored. Hence, if a BEC theory 
that reduces properly to BCS theory is at all relevant in SCs and SFs taken as many-fermion systems where pairing 
into bosonic CPs can occur, two-hole CPs must play an unambiguously crucial role. 
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